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Unitary processes allow for the transfer of work to and from Hamiltonian systems. However, to achieve 
non-zero power for the practical extraction of work, these processes must be performed within a finite-time, 
which inevitably induces excitations in the system. We show that depending on the time-scale of the process 
and the physical realization of the external driving employed, the use of counterdiabatic quantum driving to 
extract more work is not always effective. We also show that by virtue of the two-time energy measurement 
definition of quantum work, the cost of counterdiabatic driving can be significantly reduced by selecting a 
restricted form of the driving Hamiltonian that depends on the outcome of the first energy measurement. Lastly, 
we introduce a measure, the exigency, that quantifies the need for an external driving to preserve quantum 
adiabaticity which does not require knowledge of the explicit form of the counterdiabatic drivings, and can thus 
always be computed. We apply our analysis to systems ranging from a two-level Landau-Zener problem to 
many-body problems, namely the quantum Ising and Lipkin-Meshkov-Glick models. 


I. INTRODUCTION 

Recent years have witnessed a surge of interest in the study 
of thermal nanomachines that are capable of converting disor¬ 
dered forms of energy, such as heat, into useful work. At such 
small scales, thermal and quantum fluctuations play a con¬ 
siderable role, and as such the work output and performance 
of an engine are characterized probabilistically by distribu¬ 
tion functions. These distributions obey fluctuation theorems 
such as the Jarzynski equality and the Crooks equation [1-4], 
which have been verified by experiments both at the classi¬ 
cal [5-17] and quantum level [18, 19]. There have also been 
recent implementations of miniature classical engines [20-24] 
and several proposals for the realization of quantum heat en¬ 
gines [25-28]. Effects of quantum statistics of the working 
fluid have also been investigated [29-31]. 

On the other hand, it has been shown for both classical and 
quantum systems that external drivings can allow a system 
to evolve adiabatically even when driven in finite time [32- 
44]. This has applications in quantum control, and can be 
performed in three ways: (i) Driving of a system such that, 
instantaneously, a state evolves adiabatically (e.g. [32-39]). 
This is known as counterdiabatic driving and as transitionless 
driving; (ii) Protocols for which only at the final time are the 
states adiabatically transferred, while there may be excitations 
at intermediate times of the process (e.g. [40, 41]); (iii) Ap¬ 
plication of imperfect external drivings which do not allow for 
an exact adiabatic transfer, but close enough for most practical 
purposes (e.g. [42-44]). 

Hence it was suggested [39] to use such external driv¬ 
ings to render the unitary processes of a thermodynamic cy¬ 
cle quantum adiabatic while being performed in finite time. 
This could considerably augment the performance of nano¬ 
thermodynamic engines as work exchanges are extremized by 
adiabatic protocols [45]. However, implementing additional 
external driving requires resources which affect the overall 
performance of the system [46]. 

In this paper we analyze the implications of considering the 
necessary power in applying counterdiabatic driving both in a 


prototypical system such as the Landau-Zener model [47-50] 
and also in many-body quantum systems such as the trans¬ 
verse held Ising chain [51, 52]. Subsequently, we then show 
that this cost may outweigh the possible gains in work extrac¬ 
tion for slow enough processes due to the relative degree of 
adiabaticity in the dynamics. Conversely, for relatively faster 
processes, the use of counterdiabatic driving can improve the 
work exchange, depending on the experimental realization of 
the fields. Furthermore, we devise a general strategy that ex¬ 
ploits the definition of work as a two-time measurement of 
energy [53] to improve the performance of work transfer. In 
particular, we show that it is possible to achieve sizable en¬ 
ergy savings by gathering information from the first measure¬ 
ment and then applying a specifically tailored driving to the 
protocol. Lastly, we introduce an alternative measure, the exi¬ 
gency, to quantify the need for applying couterdiabatic driving 
which is related to the non-commutativity between the time 
derivatives of the Hamiltonian and the state. This measure 
has the advantage in that it can always be computed regardless 
of whether the protocol for counterdiabatic driving is known. 
Moreover, it mimics the behavior of the cost functions associ¬ 
ated to the transitionless form of counterdiabatic driving and 
goes to zero when no external driving is needed. We apply 
this measure to the analysis of the quantum harmonic oscilla¬ 
tor and the Lipkin-Meshkov-Glick infinite range spin model. 


II. COUNTERDIABATIC DRIVING 

A. The Transitionless Protocol 

To begin with, we consider the evolution of a density matrix 
p(t) from to to fi under unitary dynamics of the Hamiltonian 
Ho{t). Counterdiabatic quantum driving is obtained by apply¬ 
ing an external Hamiltonian 
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where Pj - is the projection operator on the instanta¬ 
neous energy eigenstate \j} of Hq - YjJ Ej{t)\i){i\ [32, 33, 54- 
56], 

Note that the driving in Eq.(l) is such that all energy eigen¬ 
states evolve transitionlessly. Such a strong requirement is 
however not necessary when using counterdiabatic driving to 
enhance the work output. The mean work exchange of a single 
unitary process is defined by a two time measurement proto¬ 
col = Y,m,n - En(to)] ’P’"’"pM where E,„ is 

the instantaneous energy of the m-th level, Pnih) is the prob¬ 
ability of occupying that level at to, and P'”’" is the transition 
probability from level n at time to to level m at time fi [53]. 

Next, as consequence of a resolved energy measurement 
performed before the unitary process, the system may collapse 
to one of its instantaneous energy eigenstates. Hence, given 
a particular outcome of the first energy measurement, we can 
apply a suitable driving which preserves the transitionless evo¬ 
lution of only the measured state (or a relevant sub-manifold 
of the entire system) without the need of avoiding transitions 
between the other levels that remain unpopulated throughout. 
Such a selected counterdiabatic driving for a particular eigen¬ 
state j is given by [32] 


= i^ 


dPi , 
— ,Pi 
dt 
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Here we point out that while the present discussion is pri¬ 
marily in context of unitary evolution, which is important in 
its own right, it is also relevant to the more general context 
of quantum engine cycles that consists of different strokes. 
For instance an Otto cycle is composed of two unitary strokes 
and two strokes in which only heat is transferred with the 
environment. An initial density operator that describes the 
system, under the (sequential) repetitive application of these 
strokes, generally reaches an asymptotic dynamic behavior 
which can be used to characterize the cycle. However, char¬ 
acterization of the performance of such an engine cycle in¬ 
cludes determining its net work done and/or efficiency. One 
approach to achieve this is to perform energy measurements 
after each stroke (however, we remark that in principle one 
could also determine the work done using indirect probes [57- 
59]). These measurements would then in turn affect the dy¬ 
namics of the system changing the asymptotic behavior of the 
dynamical evolution. It is thus important to also consider the 
energy measurements already within the cycle. See [60, 61] 
for a more in-depth discussion of measurement within the 
strokes of an engine cycle. 


B. Power in Generating Counterdiabatic Fields 

To quantify the power required to generate such external 
driving it is instructive to consider two examples: 

i. A single spin in a time-dependent magnetic field Bo{t), 
with Hamiltonian Ho,s - jBoit) ■ S. Here, y is the gyro- 
magnetic ratio and S, for a spin-1/2 system, is given by 
S - (/i/2) & where & is the vector composed of the Pauli 
matrices. 


ii. A neutral atom in a time-dependent electric field E(t) such 
that it experiences a potential V{t) oc |£'(f)p. 

In (i) the part of the Hamiltonian attributed to the applied 
field is given by = Bi ■ S where Bi = x (^) [33] 
which can be generated by an electric current /(f) such that 
the power required would scale as \Bi^ and thus be propor¬ 
tional to IIHfcip. However for (ii), this term could instead be 
proportional to the modulus square of another electric field 
Vzit) oc |£''(f)|^, and thus the power needed would then scale 
with the norm of the driving but not the square of it as in the 
previous case. 

Thus, while the power required to generate the counter¬ 
diabatic drivings scales as the norm of the driving Hamilto¬ 
nian (we use the Frobenius norm ||A|| = [A^Aj where A 
is an operator), the exact functional dependence on power is 
strongly affected by the experimental realization [62]. Hence, 
the cost of counterdiabatic driving can be written in general as 

c/ = Vt,n f ' \\H,rdt, (3) 

Jto 

where v,,„ is a set-up dependent constant and the index of the 
norm n depends on the nature of the applied fields [62, 63]. 
The principle of its usage here as a measure of cost is simi¬ 
lar to the constraints used in optimal quantum control studies 
[64, 65]. A closely related measure has also recently been 
used in the context of energetic cost of superadiabatic compu¬ 
tations [66]. It follows that we can define the cost of apply¬ 
ing the selected counterdiabatic Hamiltonian Nw in a similar 
fashion: 

rw.jirdt, (4) 

j 

which can be interpreted as the weighted average of the cost 
of driving each level j, over the level occupation probability 
Pj - txipPj). In fact the frequency of use of a particular driv¬ 
ing depends on the probability of measuring that particular 
energy level. Similar to v,,„, is a parameter dependent 
on the particular experimental set-up, and in the following 
we set v,„ = yw,i,n - 1 for simplicity. Here we stress that 
while (3) and (4) are entirely general expressions applicable 
to counterdiabatic fields, the exact functional dependence (n) 
and energy scale v,„ or are dependent on the particular 
form of the driving and physical nature of the fields which, 
as demonstrated by the two examples considered, cannot be 
generalized. 

We are now equipped to analyze the implications of consid¬ 
ering the cost of counterdiabatic driving in the performance of 
a work protocol. In the following, we focus on the Fandau- 
Zener model followed by the Ising model for which the coun¬ 
terdiabatic driving term is known analytically and is shown to 
be closely related to that of the Fandau-Zener model [42]. 


III. LANDAU-ZENER MODEL 

A prototypical model for the study of quantum dynamics 
is the Fandau-Zener model [47-50] which in some cases, can 
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even allow for analytical insights into the dynamics. It con¬ 
sists of a two-level system with the Hamiltonian: 

Hlz = + A(T'' (5) 


The Landau-Zener model describes the dynamics due to a 
time-dependent modulation of g(t) through an avoided level 
crossing with a finite energy gap A. For driving protocols that 
are not infinitely slow (i.e the time derivative g{t) + 0), the 
transition probability between the two levels becomes non¬ 
zero and the application of a counterdiabatic field is required 
to recover transitionless dynamics [33, 42]. Hence Eq. (1) for 
the Landau-Zener model (5) becomes 


HLZ,t - 


%(0A 

2(A2+g2)'^' 
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Note that it being a two-level system, the driving in Eq. (1) 
Hlz,! is identical to Hlz,w from Eq. (2) [67]. We now compare 
the work done in absence of any additional external driving 
(VF), with that done for an adiabatic process {Wad) considered 
in conjunction with the cost of counterdiabatic driving. The 
counterdiabatic driving is beneficial only when the cost of pro¬ 
ducing it is lesser than the inner friction (W/„e) = {W)-{Wad) 
generated in absence of the counterdiabatic fields which in¬ 
cidentally, is also the amount of additional work extractable 
attributed to the driving. 

Eor a Landau-Zener process where the state is initially far 
enough from the avoided crossing, the probability of popu- 
laton transfers decays exponentially with the timescale. This 
implies that the inner friction {Wfnc) exp[-a(fi - fo)] also 
decays exponentially, where a is time-independent. As such, 
the cost of the driving would need to decay at least exponen¬ 
tially fast with the increase of the time scale fi - fo, or the cost 
of driving would at some point be greater than the gains ob¬ 
tained by a perfect adiabatic evolution. However, we find that 
the cost of driving always decays as a power law: 


cr 


1 


(t\ - fo)" ^ 


(7) 


This can be shown by a simple change of variables s - 
(f - fo)/(fi - fo) that yields Cf = [l/(fi - fo)"-'] J ||H,(s)|rufs 
(C^ behaves analogously). It should be noted that this scal¬ 
ing is completely independent of the system and protocol 
used. Hence, the usefulness of counterdiabatic driving will 
ultimately always depend on the particular process in ques¬ 
tion. 

Despite this inherently system specific nature of the relative 
behaviors of cost and inner friction, we introduce an entirely 
general strategy that reduces the cost of counterdiabatic driv¬ 
ing regardless of the particular experimental realization. This 
involves noticing that while for a single two level system H, is 
identical to and would thus cost the same amount, they 
are vastly different for larger systems, and the resulting differ¬ 
ence in their cost can be very significant. 



FIG. 1. (Color online) (a,b) Instantaneous cost of counterdiabatic 
driving for two spins (a) and 8 spins (b). The red circles represent 
dtC] while blue continuous lines depict for decreasing values 
of the inverse temperature (as indicated by the arrow) /} = oo (top) to 
y6 = 0 (bottom). The inset of (a) shows the instantaneous exigency 
5,Co also from fi = oo (top) to y6 = 0 (bottom). 


IV. ISING MODEL 

To illustrate this, we focus on the transverse field Ising 
model because it allows for greater analytical insights with 
the exact form of the counterdiabatic field known, while at the 
same time presenting a phase transition. Its Hamiltonian reads 

(8) 

i 

In Eq. (8), J is the amplitude of spin excitation tunneling while 
g{t) is a time-dependent transverse magnetic field. In the fol¬ 
lowing we use a smooth ramp 

git) = go + (gi - go){l - cos[7r(f - fo)/(fi - fo)])/2 (9) 

in order to avoid sudden quenches [40]. 

Now, we consider the case of two spins revealing the basic 
principles of our strategy. The Hamiltonian, which we refer 
to as H 2 s, is divided into two blocks, one which dynamically 
couples the state with two spin-ups |tT) with |j,J,), and the 
other which couples Iti) with lit)- The first block reverts to 
a two-level Landau-Zener problem, which, as in (6), requires 
the application of the counterdiabatic Hamiltonian 

= "TT^T^ {ITTXUI - lUXTTI) ■ (10) 

The second block is time-independent and hence requires no 
external driving. Considering a thermal state as the initial 
condition and a change of g from 0.57 to 1.57 we observe 
a significant difference in C" and (Eig.l). In fact with a 
probability that is dependent on temperature, the first energy 
measurement could pick a state in the undriven sector which 
would thus require no external driving such that < C". In 
particular can be as small as C"/2 for small (3. 

When considering a larger system the many-body physics 
makes this study even richer. The Ising model in Eq. (8) ex¬ 
hibits a quantum phase transition at g = +7 in the thermody¬ 
namic limit, with the groundstate being paramagnetic for |g| > 
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FIG. 2. (Color online) (a) Ratio of costs Cj^/C* versus inverse tem¬ 
perature p. The arrows indicate the asymptotic value computed for 
the groundstate (at p = oo). (b) Co/C,' versus p. The blue continuous 
line represents L = 4, green dashed line L = 6 and red dotted line 
L = 8. 


J and ferromagnetic for |g| < J. Here g(t) is varied in time in 
order to cross |g| = J, without considering the case Ai —> oo 
because no counterdiabatic driving can be exactly done in that 
regime [42]. For L spins (which we take to be even in the fol¬ 
lowing), we use the Jordan-Wigner transformations [68-70] 
(jf = 1 - 2c[c,- and (xf = - (c -H cj) ]!;■<; (l - 2c^Cj) on 
to write Eq. (8) in the free-fermion form in momentum space 


Ha ^Pe 




\keke 


Pp + Pn 


\kek„ 


Po, 


( 11 ) 


where Hk - 'F/tl {[g(f) - 7cos(k)] cP - 7sin(k)cr-*} 'Fi, H(a - 
2g(l -co^co-c„^c,:). 'Fi'' = (c*l, c-k) and Pp [Po) projects on 
the even (odd) sector corresponding to the space containing 
an even (odd) number of fermionic excitations. In the even 
sector, due to anti-periodic boundary conditions, k takes the 
values ke - ±(2j - l)n/L for j e [l,L/2] while in the odd 
sector, with periodic boundary conditions, kg - +2jnjL for 
i 6 [l,L/2) in addition to = 0 and n [69, 70]. As shown 
in [42], where the authors are concerned with the groundstate, 
and focused on only the even sector, Hi k is a sum of indepen¬ 
dent Landau-Zener transitions and the counterdiabatic driving 
Hi t is thus given by 

^ /(k,f)'I';'cr^'l'i, (12) 

0<A:<;r 


where f{k,t) - -1ig{t)J sm(k)I [2^g^ + - 2g7cos(k))j and 

= {c\,C-k)- Again, depends on the state(s) selected 
upon the first energy measurement. Thus, in a similar fashion 
to the case of the two-spins, (albeit in a far richer way) the first 
energy measurement may select a reduced portion of the even 
or odd sectors; = P^^Hi tP^^, where P^^ projects over 
only the relevant set of quasi-momenta k (e.g. P^k} projects 
only over the set of even eigenstates), that are connected by 
pair creation or destruction operators (see Appendix A). As 


a result can be considerably reduced especially at 

high temperatures [see Fig. 2 (a)] [71]. The reduction in the 
energy cost of the counterdiabatic driving can be further re¬ 
duced for longer spin chains and is more significant when the 
initial condition does not include states whose dynamics en¬ 
tails small avoided crossings (in this case the cost could be 
negligible). 

Our results show that it becomes more costly to drive the 
Ising model transitionlessly as /I ^ oo (i.e. at zero temper¬ 
ature). This might seem counter-intuitive since the system 
increasingly approaches its groundstate and thus would only 
require a single eigenstate to be driven. However, while in 
general there can be considerable energy savings when the 
first energy measurement selects only a single state, the cost 
of guaranteeing adiabatic dynamics may become considerably 
large when this particular state happens to be the groundstate 
and that the system is going across a quantum phase transi¬ 
tion. In fact, it would cost an infinite amount of energy (for an 
infinitely large system) to evolve in a perfectly adiabatic man¬ 
ner [42,43]. As we consider finite systems, the corresponding 
cost for driving the groundstate will not be infinite, but instead 
be typically larger. 


V. EXIGENCY 


While C" and CJj, indicate the energy cost of achieving 
quantum adiabatic dynamics with external drivings, it cannot 
always be computed because it requires the knowledge of the 
counterdiabatic field. Moreover the counterdiabatic driving 
always ensures absence of transitions even when they would 
be perfectly balanced. For example the identity matrix com¬ 
mutes with any Hamiltonian hence it does not require any 
driving to preserve it. We thus look for an indicator for the 
need of using counterdiabatic driving which (i) would be non¬ 
zero when a driving is needed and (ii) can always be com¬ 
puted. We thus study the origin of the need of counterdiabatic 
driving by analyzing the evolution of a density matrix which is 
given by p(fi) = Up(to)U'’^ with U - T exp[-i C' HQ{t)dtl1i\, 
where T stands for the time ordering operator. Since the ini¬ 
tial conditions considered are diagonal in the basis of the ini¬ 
tial Hamiltonian Ho{to) (this includes all thermal states), to the 
lowest order in dt the evolution of p is given by 

Pit + dt) ^pit) - ^ ^ [77^^(f),p(f)] + Oidr\ (13) 

where PpQ\t) is the n-th derivative of Ho(f) and n is the low¬ 
est natural number for which the commutator in Eq. (13) is 
non-zero, details can be found in Appendix B. This implies 
that the first term of Eq. (13) which may contribute is at least 
of second order in dt. It follows that the instantaneous power 
dissipated by internal friction 'kV - fa)! dt - 0 

[72-78], where (dWfric) - {6W) - {6Wad) and (dW) is the ac¬ 
tual infinitesimal work while (dWad) is the infinitesimal work 
done if the process was quantum adiabatic. This however does 
not mean that applying an external driving to make the evolu¬ 
tion quantum adiabatic requires no power. The fact that the 
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FIG. 3. (Color online) Comparison of Exigency to the cost of transi¬ 
tionless driving in the Landau-Zener model [Eq.(5)] as g(f) is driven 
from g(to) = -lOA to g(f|) = 5A following also Eq.(9). Red con¬ 
tinuous line represents 5,C,' while the other lines show 5,Co for an 
initial groundstate occupation of pj = 1 (blue dot-dashed line) and 
Pg = 0.75 (light blue dashed line). Here t = (t - ?o)/(fi - ?o) and 

(ti - to) = A^'%. 

instantaneous variation of inner friction over time {6Wfric)l6t 
is negligible for continuous drivings implies that it is not the 
best measure of neither the need of counterdiabatic driving 
nor of its cost. Note that this is different from Ref. [79] where 
the density matrix instantaneously does not commute with the 
Hamiltonian. 

From Eq. (13), we measure the need for counterdiabatic 
driving using the quantity Co, which we refer to as exigency. 


r 

\dmt) 

L 
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Eq. (14) provides a qualitative understanding of the cost of 
driving, making it remarkably useful as it can be readily cal¬ 
culated for any Hamiltonian. Cq measure the degree of non¬ 
commutativity between p and the Hamiltonian which, if non¬ 
zero, implies the need of counterdiabatic driving. Moreover, 
(i) Co —> 0 as jS —> 0 as desired, since p is proportional to the 
identity and commutes with any time-dependent Hamiltonian 
and thus not need any counterdiabatic driving, see inset of 
Eig.l (a) and Eig.2 (b), and (ii) the instantaneous cost 5,Co, as 
shown in the inset of Eig.l (a) and in Fig.3, mimics that of the 
counterdiabatic driving 5,C/ and similarly for larger n in C". 
In particular in Eig. 3 we plot the instantaneous cost of the ex¬ 
ternal driving dtC" for a Landau-Zener problem described by 
Hiz with the same g{t) of Eig. 1 and compare it to the instan¬ 
taneous power estimated from Eq. (14), i.e. (9,Co. While these 
different measures cannot be exactly compared to each other 
because of the different, experimentally determined, constant 
v^„, they behave similarly. This means that, even without 
knowing the exact form of the counterdiabatic driving term, 
H,, it is still possible to have a qualitative understanding of 
the cost of the driving by studying Cq. We should also note, 
however, that the maximum instantaneous cost is not always 
maximum at a minimum of distance between energy levels in 
avoided crossings (as in the standard Landau-Zener problem) 
because it also depends on the exact time-dependence of the 
Hamiltonian parameters. Eor more asymmetric cases the in¬ 
stantaneous exigency <9,Co can differ, even qualitatively, from 


the cost function d,C". To further illustrate the advantage of 
the exigency we conclude with some additional examples in 
the next subsections: the harmonic oscillator and the Lipkin- 
Meshkov-Glick model. 


A. Exigency: Harmonic Oscillator 

The Hamiltonian for the quantum harmonic oscillator with 
time dependent frequency a){t) is given by 

+ f (15) 

where m is the mass and p is the momentum operator. In 
Ref [80] the exact counterdiabatic term was found to take the 
simple form 

H'^° ^--^(xp + px). (16) 

Aoj 

where ‘' ’ refers to the time derivative. Since is un¬ 
bounded C] is not finite. However the exigency can still be 
used to establish the need to perform the counterdiabatic driv¬ 
ing. In particular, we find the instantaneous power of the har¬ 
monic oscillator, as estimated by Eq. (14) for a given eigen¬ 
state reduces to 

5,Co = - {tlj\x^\il/f‘. (17) 

Eor the simple but indicative case of the groundstate this 
expression can be easily evaluated giving, \}/{x) - (x| i/r) = 

(ff)‘^'^exp(-!^^), and find {^\x^\il/) = jo?) and 
{il/\x^\4i) - Hl(2ma)). Substituting these expressions into (17) 
we finally arrive at 

5,Co = fi\d>\. (18) 

If we employ a similar ramp to that used previously, ai{t) = 
Wo + (1 - cos we can integrate Eq. (18) for f 6 

[fo, fi] and find 

Co = ^ (wi - Wo), (19) 

where the simplicity of the results stems from the choice of 
the state and the particular protocol. 


B. Exigency: Lipkin-Meshkov-Glick Model 

We now move to examining another critical many-body 
spin system, the Lipkin-Meshkov-Glick (LMG) model. A 
particularly interesting aspect is that it has infinite range in¬ 
teractions. It is therefore complementary to the short-range 
nearest neighbor Ising model previously studied. The LMG 
model can be solved analytically using the Holstein-Primakoff 
transformation [81]. The exact form of the counterdiabatic 
term in the thermodynamic limit was calculated in Ref. [43], 
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FIG. 4. (Color online) First derivative of the exigency, Eq. (14), when 
driving the groundstate of the LMG model through its quantum phase 
transition. The curves are for increasingly large system sizes from the 
bottom: = 100 (bottom, purple continuous curve), N = 200 (red- 

dashed curve), N = 300 (green dot-dashed curve) and N = 400 (blue- 
dotted curve). The top-most black thin continuous line is derived 
from the HP mapping with N = 400. Inset: The second derivative 
of the exigency. The pronounced dip that approaches the critical 
point as system size increases is clearly visible. In both hgures t = 
(t - to)/(li - to) and (f, - Iq) = 


shown to be non-local, and therefore the complexity associ¬ 
ated with engineering exact counterdiabatic driving terms was 
linked to the closing energy gap near criticality. A means to 
circumvent the requirement to implement the full correction 
term was proposed, and it was found that signihcantly less re¬ 
sources were required to achieve effective adiabatic dynamics 
when far from criticality, while more rehned correction terms 
were needed approaching the critical point. In what follows, 
we show through the use of the exigency this behavior in a 
more rigorous and quantitative manner. 

The LMG model in terms of collective spin operators, Sa - 
(where (Tq. are the usual Pauli operators) takes the 

form 


^LMG(f) = ~ (5? + yS^) - 2g{t)S,. (20) 

For the time-independent case, the LMG model has a second 
order quantum phase transition when g = A. To calculate the 
exigency we first notice that 

5Amg = -2gS,. (21) 

which allows us to directly evaluate the derivative of the exi¬ 
gency for pure states (see Appendix C) hnding 

5,Co = 2V2|g|^Var(5,). (22) 

It is immediately clear that 5,Co = 0 when g = 0 or Var(5j) = 
0. In fact this result holds for any system where the driving 
is applied to a global field. In Fig. 4 we show the behavior of 
Eq. (22), evaluated when we drive the groundstate using the 
ramp, g(f)/A = | + 5 [l - which g = A at 


t - 0.5(fi + to). Clearly, the zero points at the start and end 
of the ramp are due to g = 0. We see that as the system size 
is increased the need to apply the counterdiabatic field grows, 
and does not appear to converge except for f larger than ^(ti + 
to). Furthermore, this need grows most significantly as we 
approach the critical point, implying that the cost associated 
with driving through the quantum phase transition diverges in 
the thermodynamic limit. 

Further insight can be found by exploiting the Holstein- 
Primakoff (HP) approximation that allows us to analytically 
treat the LMG (see Refs. [43, 81] for details). This mapping 
is exact in the thermodynamic limit, and provides an accurate 
approximation for suitably large N. Setting y = 0 for simplic¬ 
ity and dropping the explicit time dependence and dehning 
g - glA for brevity, through Eq. (22) we hnd that 5,Co is 
given by 


dtCo 


2 |g| -^(gsinh^(a) + 
2 |g| sinh(a). 


0 <g < 1 

I > 1. 


(23) 


where tanh(a) = for 0 < g < 1 , tanh(a) = for 
g > 1, and we have used the HP mapping (details provided in 
Appendix D). When g > 1 we see the exigency is independent 
of the system size, however, for 0 < g < 1 the need for a 
driving term scales with increasing N. In Pig. 4 the solid black 
line corresponds to Eq. (23). The approximation has excellent 
agreement with the numerics until we approach the critical 
point, where the mapping begins to break-down for any hnite 
value of N. 

A hnal interesting point is the behavior of the second 
derivative of the exigency. In the inset of Pig. 4 we see a 
divergence that is becoming increasingly more pronounced as 
we approach the critical point for systems tending towards the 
thermodynamic limit. We remark that this behavior is equiv¬ 
alent to that of other hgures of merit which signal the emer¬ 
gence of critical behavior. 


VI. CONCLUSIONS 

Enhancing the work output of quantum engines is key to 
designing future nano-technologies. External helds for coun¬ 
terdiabatic driving can increase the efficiency of thermal ma¬ 
chines by increasing the work extracted in the unitary strokes. 
However the effectiveness of this method strongly depends on 
the cost of applying the external driving and on the duration of 
the process. While a quantitative evaluation is dependent on 
the particular experimental realization, in general, at longer 
time-scales the use of counterdiabatic driving is detrimental 
because the energy required to generate the held is larger than 
the energy gained from a quantum adiabatic evolution. The 
cost of applying counterdiabatic driving can be signihcantly 
reduced by choosing a selected form of driving which depends 
on the hrst measurement of energy, especially in systems close 
to a phase transition where the exact external driving is par¬ 
ticularly costly for states more strongly affected by the tran¬ 
sition. In future, optimizations comparing perfect against ap- 
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proximate counterdiabatic drivings and hybrid protocols that 
selects perfect or approximate drivings depending on the ini¬ 
tial energy measurement could be implemented for consider¬ 
able energy savings especially as quantum technologies are 
scaled-up. 
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Appendix A: Case of 8 spins 


Here, we consider the illustrative case of 8 spins for the 
Ising model explicitly. In this scenario, the Hilbert space com¬ 
prises of 256 states, 128 each in the even or odd sector. In the 
even sector the possible values of the quasi-momentum are 
ke = ±7r/8, +3;7r/8, +57r/8, +77r/8, while in the odd sector 
ko = 0, +7r/4, +7r/2, +37r/4, n. 

With a little computation, it becomes apparent that both 
even and odd sectors of the Hamiltonian are further divided 
into sub-blocks. This is due to the fact that in the basis 
spanned by the operators c\ acting on their vacuum |v) (c],|v) = 
0 for every k), only states differing by a pair of creation or de¬ 
struction operators are coupled dynamically. 

Eor instance, one sub-block of the even sector, which we 
will refer to Sa, is spanned by the 16 basis elements 


Sa ' 


|v) 

p\\v),pI,\v),pI\v),pI\v) 
p\pI\v), p\pI\v), pIpI\v), pIpI\v), pIpI\v), 
p\pIpI\v}, pIpIpI\v), pIpIpI\v), pIpIpI\v) 
p\pIpIpI\v) 


where we have adopted the notation p^. - . 

“8 “T 

Another block, S b, is formed by single states that are 
completely uncoupled to any other state, (for example 
ctc^ 3 „C 5 „c^, |v)) because no pair p\ can be added to or be re- 

8 T T T ■' 

moved from it. By simple combinatorics, it is apparent that 
there are 16 such states in the even sector. Eor example 
ctc^ |y) is also one such state that can be obtained 

8 ~T ~T T 

combinatorially. 

A qualitatively intermediate scenario occurs in the sub¬ 
blocks Sc- Eor example one sub-block of 5c could contain 
the states that are spanned by c\c\^ |v), c],c\^pl\v), c\c\,pl\v), 

8 t 8 T ST 

c\c\^plpl\v}. There are 24 such sub-blocks in the even sector 

8 T ^ 

of the Hamiltonian that are of a similar structure (for exam¬ 


ple another group of 4 coupled states is given by 

“T T 

c\j,p\c\jv), p\c\j\jv), p\c\^pIc\Jv)). 

8 8 8 8 8 8 

Now putting it all together, we find that the sub-block S a 
has 16 states, all the sub-blocks of type Sb have a total of 16 
states while there are a total of 96 states in the sub-blocks of 
type S c- The sum of all these states is indeed 128 as expected. 
A similar scenario unravels in the odd sector. 

It is now clear that if the first energy measurement selects 
a state in any of the S b sub-blocks, no transitionless driv¬ 
ing is needed because the states in 5 b are invariant during 
the time evolution. If instead the measurement selects an 
eigenstate of the sub-block Sa, then, in order to keep the 
evolution transitionless, it will be necessary to apply ///„, = 
Pe Pe with k given by all the possible kg. 

Lastly, for a state in 5 c, it would be sufficient to drive only two 
values of k (which specific values of k to be driven depends on 
which pairs are involved in the sub-block). Eor instance, in the 
example above, k = 57r/8 and In/S are needed. 


Appendix B: Time Evolution of a Diagonal Density Matrix in 
the Instantaneous Basis of the Hamiltonian 


Here we furbish the details in deriving Eq. (13). Assuming 
that the series expansion of p(t + dt) converges for sufficiently 
small dt, we write 


p(f -H dt) 


y d”p(t) dt" 
^ dt" n\ 


Erom here, using 


dp(t) 

dt 


[^o( 0 .p( 0 ] 


we obtain the higher derivatives of p, for instance. 


d^p 

dfi 




cm 

dt 


(Bl) 


and 


d^p 

dP 



Ho, [Ho, 


Ho, 

dHo 


i \d^HQ „ 



(B2) 


Hence, it follows that if the n - 1 derivatives of Ho commute 
with p, then the lowest order correction in dt to the time evo¬ 
lution of p will be given by Eq. (13) in the main text. 
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Appendix C: Exigency for the LMG model 


Here we detail the calculation to arrive at Eq. (22) in the 
main text. Assuming the state we wish to drive is p = li/r) {i/f| 

dfCo - ||[p, 5,i/LMG]|| 

= \\p(-2gS,) - (-2gS,)p\\ 

- 2\g\ ^Tr[(S,p-pS,)(pS,-SS] 

= 2-\g\ -y/xr [S JP-S^-SMzP-pS-Mz+pSI p\ 

= 2V2|g|^Tr[p252-5^5,p] 

- 2 V2|g| 

= 2V2|g| Vvar(5,). 


Appendix D: Exigency for the LMG model: Holstein Primakoff 
approximation 

Here we outline the steps required in order to derive 
Eqs. (23), with j - 0 for simplicity. This requires us to cal¬ 
culate Eq. (22) using the Holstein-Primakoff (HP) transfor¬ 
mation. Eor suitably large N we map the spin operators into 
the creation and annihilation operators a and of a harmonic 
oscillator 


— [a + a'j, 

(Dl) 

N 

- a a. 

2 

(D2) 


Eollowing the Supplementary material of [43], when g > 1 the 
HP transformation is always taken along 5^. In order to map 
the LMG model to the harmonic oscillator we are required to 
perform a Bogoliubov transformation 

a = sinh(f)/7l+cosh(f)Z7, (D3) 

= sinh(f)Z7 + cosh(f)Z7'', (D4) 


with tanh a - Therefore to calculate the required expec¬ 
tation value we must write a^a in terms of b and b^, and cal¬ 
culate the expectation value over the groundstate. It is readily 
found that the only term (other than constants) contributing is 
proportional to b^b^^. As such, we find that the variance of 5^ 
is 


Var (^z) = ^ sinh^ (a). (D5) 

Substituting into Eq. (22) we arrive at Eq. (23) (b) 

5,Co = 2|g| sinh(a). (D6) 


A similar, albeit more involved, calculation is required for 
0 < g < 1. In this case the direction along which the HP 
transformation must be taken changes with the value of g. In 
this case the operator of which we must calculate the variance 
of is 


= ^jcosycj + 5xSinyf), (D7) 

with cos (f - g. Once again, we express 5 f in terms of b and 
b '^ with tanh a - , and calculate the expectation value over 

the groundstate. We now find that the only terms contributing 
are constants, and terms proportional to bb'^ and hP'b^^. After 
some manipulation we find 


Var(5j) = ^1 sinh^(a) + ^ [sinh(|) + cosh(|)j^j 

= ^1 sinh^ (a) + . 

(D8) 

Substituting into Eq. (22) we arrive at Eq. (23) (a) 

5,Co = 2|g| ^i^gsmhHa)+^^^e-j (D9) 
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